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СОМУ5ТКОСТІХС ВЕ5Т СНЕВУ5НКЕУ 5РІЛХЕ АРРКОХІМАТІОМ5 


З | метою побудови найкращого  чебишовського - наближення для заданої функції 
поліноміальним сплайном степеня п з "г фіксованими вузлами у статті пропонується застосувати після 
відповідної модифікації алгоритм апроксимації функції багатьох змінних узагальненим многочленом. 
У цьому алгоритмі використовується зведення до задачі лінійного програмування з головною 
двоїстою максимум-задачею. Аналіз чисельних результатів показав, що у більшості випадків 
модифікований алгоритм знаходить більш точні наближення сплайнами, ніж інші відомі алгоритми. 

Ключові слова: поліноміальний сплайн, найкраща чебишовська апроксимація, алгоритм, 
лінійне програмування. 
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Вступ 

Для апроксимації функціональних залежностей складної структури, які 
виникають при розв'язанні різноманітних  науково-технічних задач, широко 
використовують сплайни. Найважливішу роль серед них грають поліноміальні 
сплайни, склеєні з кусків многочленів. Вони застосовуються, наприклад, для 
аналітичного опису поверхонь деталей машин (крила та фюзеляжу літака, корпусу 
водного судна, профілю сопла реактивного двигуна), для завдання траєкторії руху 
різця станка з програмним керуванням, при розрахунках на міцність будівельних 
конструкцій та ін. 11-31. 

У статті розглядається задача побудови поліноміального сплайна найкращого 
чебишовського наближення з фіксованими вузлами. Цей сплайн при такому ж числі 
параметрів наближає функцію не гірше, ніж інтерполяційний сплайн. Крім того, для 
побудови інтерполяційного сплайна звичайно потрібно задавати ще деякі граничні 
умови. Тому в багатьох випадках сплайн найкращого чебишовського наближення 
виявляється кращим |41. 

Постановка задачі, аналіз останніх досліджень і публікацій 


Нехай на скінченому відрізку (а, В| дійсної осі задано сітку 
А, Су СУрЗ--ЧУ, ЧУ ТВ, (1) 
де г -- натуральне число. 


Функцію х(х) називають поліноміальним сплайном степеня п дефекту К 
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(1«Кк«п) з вузлами  уі,...У,, якщо на кожному відрізку ГУ ул М 
(7 -1,...,г-1) функція 5 збігається з деяким многочленом степеня п 1 
зеє""е рі: 


Прикладом сплайна степеня п дефекту І з єдиним вузлом 5 може бути функція 
(х-6)",х»6, 
0, же. 


Множину сплайнів (5(х)) для фіксованої сітки Д, позначимо 7, а, Це - 


скоб з (2) 


лінійний простір розмірності п-кук--І. Будь-який сплайн 5є7,,л МОЖНА 


представити у вигляді суми многочлена і лінійної комбінації елементарних сплайн- 
функцій (2), причому таке представлення єдине |3|. Зокрема, сплайн степеня п 
дефекту І записується у вигляді 


п є ї 
502) з У а (х-а) -Уара(х-міг» (3) 
ід ічІ 
де ад д.д, 7 Дійсні числа. 


Задача найкращого чебишовського наближення заданої на відрізку |а, ВІ 
функції /(х) поліноміальним сплайном з г фіксованими вузлами (1) полягає у 


. . ж - : 
знаходженні сплайна 5 є 7, кл, ЯКИЙ ЗадоВОЛЬНЯЄ умові 
ж 2 
тах |/(9)-5 (х) 2-2 о тіп тах ГО - 509) ; (4) 
хе, ВІ 509є7п, ЮА, хе а, В 
Труднощі розв'язання задачі (4) пов'язані, насамперед, з тим, що простір 
сплайнів 7, ;л, Не Є чебишовським. Тут не працюють напряму наближені методи 


су ж Ч - є 
побудови сплайна 5 , які грунтуються на теоремі Чебишова |б6). Крім того, у 
загальному випадку може існувати більше одного сплайна найкращого наближення. 


Для знаходження сплайна 5 у роботах |4,7, 8| пропонується замінити відрізок 
сь В| деякою сіткою 


Ап! Фехяєх є. «х,оеВоетьпяг, 
і звести задачу найкращої чебишовської сплайн-апроксимації на цій сітці 
ж . 
тах | /(х))- 5 (х)Д 2; опіпо пах | /(ху) -5(х)) (5) 
1«І«т ЧЕ УДЛИЧИ 1«І«т 


до задачі лінійного програмування (ЛП) 
пі 2, 
70) - Хай з а) -Уан(жр зУД)"-А. 40,1 21,..., т, (6) 
і-0 чі 


-гокУ або зад Уам (о зу) -А «0,121... т, 
ї«0 


ічі 
є ; РАрІА 8 з ж 

методи розв'язання якої добре розроблені. Вектор коефіцієнтів сплайна 5 (х) можна 

також шукати як чебишовську точку несумісної системи лінійних рівнянь |4). Третій 


Й і - ж ; й 
підхід до знаходження сплайна 5 (х) полягає в узагальненні алгоритму послідовних 
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чебишовських інтерполяцій Ремеза |9--11| на випадок апроксимації поліноміальними 
сплайнами (12, 131. 

У статті пропонується розв'язувати задачу сплайн-наближення (4) зведенням 
до алгоритму А.О. Каленчук-Порханової та Л.П. Вакал (14) побудови узагальненого 


М 4 М М 
многочлена Ф (х) - Уа; Ф, (Х ) найкращого чебишовського наближення для функції 
і-0 


ЛОЮ), Хей(х,оьх,)є С, де О - множина точок р-вимірного евклідового 
простору, (Ф.-(Х)|4 о - система лінійно незалежних на О функцій, які називаються 


7 ж - . 
базисними. Узагальнений многочлен Ф (х) є розв'язком наступної задачі 


пах ЦСК) -Ф'ЄХ,) - піп тах СХ) - У7амфК. б 
Хео а (Хе0О 20 


В алгоритмі (14| використовується зведення задачі (7) до задачі ЛП, при цьому, 
завдяки урахуванню специфіки задачі чебишовської апроксимації, вказаний 
алгоритм є більш ефективним для її розв'язання, ніж загальні методи ЛП (див. 
докладніше в наступному розділі), У випадку апроксимації поліноміальним 
сплайном виду (3) як базисного в алгоритмі вибирають функції 

ФОд 2(х-а)", і -1,...п-Ї; Ф.О2 2(х-м,,аї інтяє?,патяі. 

У низці прикладних задач зручно використовувати представлення сплайна 

5(х) через так звані В-сплайни, спеціальним чином побудовані функції М,,(х) з 


2 кл, » ЯКі ВИЗНачають за рекурентними формулами (21: 


1, хе |М, Уа Х-Мі 


Мода М мл УНН 9 8) 
0, х ГУ Угці і м, -Мм й і 


іп і Уйнтчі 7 Уїні 


У од з 


Функції М, (х) (і є -пу-п -1,...,г ) є лінійно незалежними. Якщо ввести додаткові 


вуши У «У зно чУДхО5 Вона Ура без то будь-який сплайн 


гепя22 
5 Є 27 ,ха, МОЖНА оДНОЗНачно представити у вигляді |6| 


пат 


коре ЗаМозьїКої: (9) 
сі 


Алгоритм (14| апроксимації узагальненим многочленом після відповідної 
модифікації може бути застосований і для знаходження найкращого чебишовського 
наближення функції /(х) на множині точок сітки Д,, за допомогою сплайна (9). У 
цьому випадку в ролі базисних виступають функції (8). 

Модифікація алгоритму для апроксимації сплайнами 

В алгоритмі |14| реалізовано підхід, в якому розв'язання задачі найкращого 
чебишовського наближення (7) для функції /(Х), Х нор на деякій р- 


з 4 
вимірній сітці узагальненим многочленом Ф (х) - Уа Ф; (Х ) за системою базисних 
ізо 


функцій (Ф.(Х))7 », зводиться до задачі ЛП з головною двоїстою максимум-задачею, 
число обмежень якої при т»» 4 (т- число точок сітки) значно менше кількості 


обмежень прямої задачі. У випадку застосування алгоритму для апроксимації 
сплайнами виду (3) або (9) число обмежень двоїстої задачі дорівнює п--г-2,а 
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прямої - 2гт-(п-гч2). 

Алгоритм |14| дозволяє отримати найкращі наближення з меншою похибкою 
та, як правило, за значно меншу кількість кроків, порівняно з аналогічними 
алгоритмами (наприклад, алгоритм |15|). Це досягається за рахунок низки 
удосконалень обчислювальної схеми алгоритму, основним з яких є такі: 

- зведення до максимум-задачі ЛП проводиться з орієнтацією на отримання 
початкової симплекс-таблиці з якомога більшим значенням цільової функції, 
завдяки чому оптимальний розв'язок знаходиться за відносно малу кількість кроків; 

- стандартна  симплекс-таблиця замінюється майже удвічі меншою, але 
рівноцінною за поданою інформацією таблицею, в якій вже є допустимий 
базисний розв'язок двоїстої максимум-задачі; 

- при розв'язанні максимум-задачі ЛП модифікованим  симплекс-методом 
реалізується напівоптимальний варіант переходу від одного допустимого 
базисного розв'язку до наступного. 

При модифікації алгоритму (14) для апроксимації поліноміальними сплайнами 
потрібно змінити набір вхідних даних. У цьому випадку для роботи алгоритму слід 
ввести такі дані: 

- кількість коефіцієнтів сплайна 4 (д-п-г-і); 


- вузли сплайна у,, / - 0,...,Г-1 (мо 2 4,У, -В); 
- число точок сітки т; 

- точки х; сітки, І - 1,..., т; 

- значення функції у точках сітки /(х,). 


Вважається, що т» п--г--1 і на кожному з проміжків |у ізУЇ, )21,...г41, 


є щонайменше дві точки сітки Д,,. 


Для побудови сплайна потрібно також внести зміни в процедуру формування 
матриці, на базі якої створюється опорна симплекс-таблиця. Ця матриця 
позначається Т 1 має розмір (п- г -3)хт. У випадку апроксимації сплайном виду 


(3) елементи 1; ; цієї матриці обчислюються за формулами: 
іо зігдє 
іц - (хз - 0), Р21,,..п41, 152 (0 Ура» Ї ПА 2,аПАГУІЇ, 


їла, і з -Го;) , Їлгя3,) -0,)-1,..,т, 
Якщо ж для наближення використовується представлення сплайна у формі (9), 
то елементи матриці Т визначаються за формулами 
р ЗУ (0) хо, ЗЛО) ної 051, 2-12, 
Для обчислення значення В-сплайна М,,(х) степеня п в точці х з проміжку 
ГУ г» їм | розроблена спеціальна процедура Вуріїпе, розрахунки в якій проводяться з 
використанням рекурентних формул (8). При цьому додаткові вузли, необхідні для 
обчислення В-сплайнів, покладаються рівними а, та В, тобто му , -..-уУ |за, 


Ура ПЗ Уачо З В. Число операцій множення і ділення, що виконуються в 


процедурі Вуріїпе для кожного значення х, пропорційне п! 
Результатами роботи модифікованого алгоритму є значення похибки р 


ТІ зб - є ж 
найкращого чебишовського наближення заданої функції / сплайном 5 
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р епих | /(х/)-5 (ху) 


1«І«т 


Аа Х М ж ХУ 
та значення коефіцієнтів ад й» ..,а пр сплайна 5 , якщо використовується сплайн У 


пз 


ц» ЯКщО СпЛайн представляється у формі (9). 


формі (3), або коефіцієнтів сі,....Сууч 

Результати обчислювальних експериментів 

Для перевірки ефективності застосування модифікованого алгоритму для 
апроксимації сплайнами з фіксованими вузлами виконано серію обчислювальних 
експериментів. Для низки тестових функцій, результати сплайн-апроксимації яких 
відомі з наукової літератури, похибки наближення, обчислені за цим алгоритмом, 
порівнювались з похибками, знайденими за іншими алгоритмами найкращого 
чебишовського наближення сплайнами. 

У табл. 1 1 2 представлені результати для двох тестових прикладів, а саме: у 


табл. 1 наведено похибки апроксимації функції /(х)-х на проміжку (0,1) 
сплайнами степенів від другого до п'ятого, а в табл. 2 - похибки наближень функції 
лозаз ху" на відрізку 10, || сплайнами третього 1 четвертого степенів, отримані 
за модифікованим алгоритмом (у таблиці позначається як алгоритм 1) і за 
алгоритмом (13) (позначається як алгоритм 2). Зазначимо, що для обчислень у 
першому прикладі використовувалась рівномірна сітка з т - 2001 точок, у другому - 
зт 1001 точок. 


Таблиця 1. Апроксимація функції /(х) - Ух на Г0, П) сплайном степеня п 


Число Вузли Похибка апроксимації 

вузлів г па2 па3 па-4 паб 
1 0,04 0,04077 0,01228 0,01041 0,00917 
2 0,0065; 0,108 0,01190 0,00495 0,00402 0,00332 
3 0,002; 0,02; 0,15 0,00840 0,00312 0,00181 0,00137 


Таблиця 2. Апроксимація функції /(х) (1 х)"" на Ї0, П сплайном степеня п 


Число Вузли по 3 па4 
вузлів Похибка за Похибка за Похибка за Похибка за 
г алгоритмом 1 | алгоритмом 2 | алгоритмом І | алгоритмом 2 
1 0,5 0,000325 0,000385 0,000051 0,000051 
2 0,33; 0,67 0,000086 0,000089 0,0000099 0,000010 
3 0,25; 0,5; 0,75 0,000033 0,000033 00000036 00000036 
4 0,2; 0,4; 0,6; 0,8 0,000015 0,000015 0,0000017 0,0000017 
Висновки 
Для побудови  поліноміального - сплайну найкращого  чебишовського 


наближення, у статті запропоновано застосовувати, після відповідної модифікації, 
алгоритм апроксимації функції багатьох змінних узагальненим поліномом |14). В 
алгоритмі використовується зведення задачі чебишовського наближення до задачі 
лінійного програмування з головною двоїстою максимум-задачею, число обмежень 
якої значно менше кількості обмежень прямої задачі. Завдяки удосконаленню 
обчислювальної схеми цей алгоритм має переваги за точністю та швидкодією у 
порівнянні з іншими аналогічними алгоритмами. 

При модифікації алгоритму для наближення сплайнами розширено набір 
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вхідних даних, змінено процедуру формування матриці, на базі якої створюється 
опорна симплекс-таблиця, додано процедуру обчислення значень В-сплайнів у 
точках (для випадку подання поліноміального сплайна у вигляді лінійної комбінації 
В-сплайнів) та їн. Аналіз результатів обчислювальних експериментів показав, що в 
більшості випадків за допомогою модифікованого алгоритму вдається побудувати 
сплайн-наближення з кращою точністю, ніж за іншими відомими алгоритмами. 
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КЕ5ОМЕ 


І, Р. Уакаї 

Соп5ігисбіпе Безі спебузПпеу 5ріїпе арргохітайоп5 

Бог а біуеп Гипсйоп /є а, | ууе аге іп(егезісд іп соптрийпеє (Фе Безі 
арргохітацоп їтот Ше 5иберазе ої зріїпе Гипспопз5 ої дергее п апа деїесі К «І міф г 
Нхеад Кпоїз а «му «...«м, «В, УЮ гевресі го Ше Сребузреу (5иргетит) погт. 

Іо огдег іо соп5ігисі (пе Бе5і зріпе арргохітапоп, 1ї( 15 ргороз5еай (о арріу, айгег ап 
арргоргіаїе тодійсацоп, ап аїбогійт їог арргохітайтя плапу-уагіабіе5 Гипсйоп Бу а 
бепегайгед роїупопаиа! (14). Гп Фе аїсогійдті а гедисйоп (о Ше Шпеаг ргоггаттіпе 
ргобіета ул Фе паїп иа! плахітита-ргобіет 15 ц5ей. ТБе аїєогійрт 15 плоге еНесПує 
Фап єепега! тешодз ої Ппеаг ргобгатитіпя 5їпсе 18 (акез їпіо ассойпі Ше рагіїсиіагігу ої 
Фе Без: Спебу5реу арргохітайоп ргобіет. 

Тріз аїєогійт їог арргохітайтпеє папу-уагіабіез Гапсйоп у/аз5 плодійеа їог зрійпе- 
арргохітацоп. і рагіїсиіаг, Ше іприї 5е( у/а5 ехрапаед, Ше ргоседиге ої таїйтіх 
Рогтайоп (оп Фе Ба5і5 ої Ші8 тайтіх а 5ітрієх їабіе 15 сгеагед) має срапеедй, а ргоседйге 
Гог согарийпе В-5ріїпез уаїше5 м/аз5 аддед (мреп а 5ріїпе 15 ргезепісд Бу а Ппеаг 
сопабіпайоп ої В-5ріїпев5), еїс. 

Тре тодібеа аїсогіїт Па5 Бееп (езіса оп а питабег ої Гапсййоп5. ГП Фе агіїсіеє ме 


яїме бо ої Ше5е питегіса! ге5ціїв. ТРре бг58: іпуоЇуез арргохітайте, Ух оп ГО,П Бу 


зрійпез ої деєгее5 п -2,5 апі Фе зесопі - (Ін х) | оп (0,1) Бу взріїпез ої дергее5 
па3,4. Апаїу5і5 ої Ше питегісаї ге5иціїє 5роууед Ша іп піо5і са5е5 Ше подіПпед 
аїбогіййт сотриїез5 зріте арргохітайопя поге ргесі5еїу Фап обег Кпом/п аїєогіПт8. 
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